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ABSTRACT 

We present Spitzer 5-15 fim spectroscopy of one cool white dwarf and 3.6 - 8 
/im photometry of 51 cool white dwarfs with T e g < 6000 K. The majority of our 
targets have accurate BVRIJHK photometry and trigonometric parallax mea- 
surements available, which enables us to perform a detailed model atmosphere 
analysis using their optical, near- and mid-infrared photometry with state-of-the- 
art model atmospheres. We demonstrate that the optical and infrared spectral 
energy distributions of cool white dwarfs are well reproduced by our grid of mod- 
els. Our best fit models are consistent with the observations within 5% in all 
filters except the IRAC 8 /im band, which has the lowest signal-to-noise ratio 
photometry. Excluding the ultracool white dwarfs, none of the stars in our sam- 
ple show significant mid-infrared flux deficits or excesses. The non-detection of 
mid-infrared excess flux around our f«2-9 Gyr old targets constrain the fraction 
of cool white dwarfs with warm debris disks to 0.8±oJ %• 

Subject headings: infrared: stars - stars: white dwarfs - stars: atmospheres 

1. Introduction 

Devoid of a ny nuclear bu rning, white dwarfs simply cool with time. This well understood 



cooling process (lMestellll952l ) enables astronomers to use them as cosmic chronometers. The 



absence of white dwarfs fainter than log(L/L ) = —4.5 in the Solar N eighborhood is used 



as evidence of the finite age of the Galactic disk (IWinget et al.l 119871 ). The best estimate 



for the age of the Galactic disk based on the white dwarf luminosity function is 8 ± 1.5 Gyr 
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( ILeggett et al.l Il998l ; lHarris et al.l 120061). This method has no w been applied to two halo 
globular clusters, M4 and NGC 6397 Jrlansen et al.lhood boOTT ). 



Prior to the Spitzer Space Telescope, the oldest white dwarfs in the Solar Neighborhood 
could not be observed in the mid-infrared due to their faintness. Therefore, the luminosi- 
ties of cool white dwarfs were estimated using optical and near-infrared photometry, and 
bolometric corrections based on white dwarf model atmospheres. An ob servational c heck o n 
these corrections is required to c onfirm their reliabil ity. Such studies by lKilic et al.l (j2006al ); 
Tremblay fc Bergeron! (120071 ) and lKilic et al.l (]2008al ) showed that the white dwarf models are 
able to explain the spectral energy distributions (SEDs) of almost all cool white dwarfs with 
^eff > 6000 K. However, significant mid-infrared flux deficits compared to model predictions 
were discovered for one star, LHS 1126. 

The main atmospheric constituents in cool white dwarfs are hydrogen and helium, 
though trace amounts of heavier elements are also observed in some white dwarfs. The 
atmospheric compositions of white dwarfs hotter than 5000 K can generally be constrained 
based on the presence of Balmer lines. However, Balmer lines disappear below about 5000 K, 
and optical spectroscopy cannot be used to identify hydrogen versus helium rich white dwarfs. 
The most critical opacity sources in cool hydrogen-rich w hite dwarf atmospheres are believed 
to be collision induced absorptio n (CIA) in the infrared ([Hansen 1 1998 



Saumon &: Jacobson 



1999h and Ly a in the ultraviolet fcowalski fc Saumon&d : iKoester fc Wolflfeoooh. The pri- 



mary opacity source in helium-rich white dwarfs is He free- free absorption (IKowalski et al. 



20071 ). Helium atmospheres are much denser than hydrogen atmospheres at the same temper- 
ature. Therefore if hydrogen is present, collisions between H2 molecules and neutral helium 
can induce CIA. This opacity becomes significant at higher temperatures compared to pure 
hydrogen atmospheres. 

CIA opacity is expected to produce molecular absorption features in the near-infrared, 
but current CIA calculations do not predict absorption bands in the mid-infrared. Therefore, 
the observed mid- infrared flux deficits for LHS 1126 imply that either the CIA calculations 
in white dwarf atmospheres are incomplete, or these flux deficits are caused by some other 
mechanism 0. An important question for white dwarf cosmochronology is whether LHS 1126 
is unique or not. The previous Spitzer observations of cool white dwarfs included only a few 
stars cooler than 6000 K, and the mid-infrared flux distributions of the oldest white dwarfs 
in the Solar Neighborhood have not been studied yet. 



In order to understand the mid-infrared SEDs of cool white dwarfs and the CIA opacity, 



1 LHS 1126 shows blue shifted carbon bands in the optical fsee lHall fc Maxwellll2008l ) and we cannot rule 
out peculiarities in its infrared SED due to the presence of carbon or other elements. 
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we used the Infrared Array Camera (IRAC, iFazio et al.ll2004j ) on the Spitzer Space Telescope 
to observe 53 nearby cool white dwarfs. We also used the Infrared Spectrograph (IRS) to 
observe one of our targets. In this paper, we present IRAC photometry and IRS spectroscopy 
of these stars and perform detailed model atmosphere analysis using both optical and in- 
frared data. Our observations and model fits are discussed in §2 and §3, respectively, while 
implications of these data are discussed in §4 and §5. 



2. Sample Selection, Observations, and Reductions 



We selected our targets from the samples of iBergeron et al.l (119971 . 1200 ll . 120051 ) . The 
beauty of this sample selection is that Bergeron et al. provided accurate BVRIJHK photom- 
etry and trigonometric parallax measurements for the majority of our targets. Our sample 
also includ e s three cool white dwarfs found in the Sloan Digital Sky Survey (SDSS) by 
Kilic et al.l (j2006bl ). The se three stars hay e near-infrared photometry from the Two Micron 
All S ky Survey (2MASS , ICutri et al.lliooih. We also add ed three ultracool whi te dwarfs, LHS 



3250 flHarris et alJll999h . CE 51 flRuiz fc Bergeronll200lh . and WD 0346+246 flHamblv et al. 



19971 ). to our sample to extend our study to T e g < 400 K. The majority of o ur targets have 
accurate coordinates and proper motions measured by iLepine fc Sharal (120051 ) , which is very 
helpful in identifying these faint objects in crowded fields. 

Observations reported here were obtained as part of our Cycle 3 GO-Program 30208 (PI: 
M. Kilic). We obtained 3.6, 4.5, 5.8, and 8 /mi images for 53 cool white dwarfs including 
two common proper motion systems in 51 Astronomical Observation Requests (AORs). 
Depending on the source brightnesses, integration times of 30 or 100 seconds per dither, 
with five or nine dithers for each target, were used. We reduced the data using both IRAF 
and IDL r o utines. Our reduction pro c edure s are similar to the procedures employed by 
Kilic et al.l J2008ah and iMullallv et all (l2007h . Briefly, we use the IRAF PHOT and IDL 



astrolib packages to perform aperture photometry on the individual BCD frames from the 
latest available IRAC pipeline reduction (S14.4.0, S15.3.0, S15.0.5, or S16.1.0 for our targets). 
In order to maximize the signal-to- noise ratio, we used a 5 pixel aperture for bright, isolated 
objects, and 2 or 3 pixel apertures for faint objects or objects in crowded fields. Following the 
IRAC calibration procedure, corrections for the location of the source in the array were taken 
into account before averaging the fluxes of each of the dithered frames at each wavelength. 
Channel 1 (3.6/mi) photometry was also corrected for the pixel-phase- dependence (see the 
IRAC Data Handbook). The results from IRAF and IDL reductions were consistent within 
the errors. The photometric error bars were estimated from the observed scatter in the 5 
(or 9) images (corresponding to the dither positions) plus the 3% absolute calibration error, 
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added in quadrature. Fin ally, we divided the estimated fluxes by the color corrections for a 



Rayleigh- Jeans spectrum (jReach et al.ll2005l ). These corrections are 1.0111, 1.0121, 1.0155, 



and 1.0337 for the 3.6, 4.5, 5.8, and 8 /zm bands, respectively^. 

We present the average fluxes measured from the Spitzer images in Table 1 . Two of the 
stars in our sample, WD 0346+246 and WD 1310—472 were blended with brighter stars, 
which prohibited us from performing photometry on these objects. These two stars are not 
included in Table 1. Another common proper motion binary system, WD0727+482A and 
WD0727+482B, was unresolved in IRAC images, and the photometry presented in Table 1 is 
the total flux from the system. We included this flux measurement for completeness purposes, 
however we do not include these two stars in our analysis for the remainder of the paper. 
Five other stars, WD 0222+648, WD 0551+468, WD 0851-246 (CE 51), WD 1247+550, 
and WD 2054—050, are near brighter stars in the IRAC images, and the photometry is 
probably affected by these nearby sources. We present the photometry for these five stars 
and include them in our analysis, but do not use them to draw any conclusions. 



As part of our Cycle 3 program, we also used the Infrared Spectrograph (IRS. lHouck et al 



20041 ) Short-Low module in the first and second orders to obtain 5-15 /im spectroscopy of 
the brightest cool white dwarf in our sample, Wolf 489 (WD 1334+039). These observations 
consisted of 44 cycles of 1 min exposures and were performed on 2007 July 29. For each of 
the 44 basic calibrated spectral images (BCDs) at each of the two nod positions and each 
of the two orders of the IRS, a cube was compiled. Each cube was collapsed into a spec- 
tral image, using a robust average. Outlier pixels with respect to a 5 x 5 box median were 
replaced with the median, to account for rogue and dead pixels. The two nodded spectral 
images for each order were differenced. A significant spatial gradient is evident in the nod 
difference images, so we removed a row-by-row (i.e. approximately constant wavelength) 
median spanning spectral order on the nod- dif ference image. The background-subtracted 



difference images were analyzed using SMART (IHigdon et al. 1120041 ). using manual extrac- 
tions of tapered apertures and local background subtraction from the portion of the slit not 
covered by the source. The extracted spectra were merged and rebinned using a logarithmic 
bin spacing, yielding the final spectrum. 



3. Model Atmospheres and the Fitting Procedures 

We use state of the art white dwarf model atmospheres to fit the BVRIJHK and Spitzer 
photometry of our targets. Our fitting procedure is the same as the method employed by 



These corrections were also used in iKilic et al.l (|2006af) and iKilic et al.l (|2008af ) . 
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Tremblay fc Bergeronl (120071 ). Briefly, the magnitudes are converted into monochromatic 
fluxes using the zero points derived from the Vega (STIS) spectrum integrated over the pass- 
band for each filter. The resulting fluxes are then compared with those predicted from the 
model atmospheres, integrated over the same bandpass. We use the parallax measurements 
to constrain the surface gravity, if available. Otherwise, we assum e a surface gravity of log g 



8. The model atmospheres include the Ly a far red wing opacity (IKowa 



as we ll as non-ideal physics o f dense helium that i ncludes refraction (IKowalski fc Saumon 



ski&Saumon 2006) 



20041 ). ioniz ation equilibrium (IKowalski et al.l 12007), and the non-ideal dissoc iation equilib 



rium of H 2 (IKowalski! 120061 ) . Using these models, IKowalski fc Saumon! (120061 ) were able to 
reproduce the ultraviolet to near-infrared spectral energy distributions of several cool white 
dwarfs including the DA star BPM 4729. This star has a UV spectrum that can only be 
reproduced by pure hydrogen atmosphere models including the Ly a opacity. This opacity 
occurs due to collisions of hydrogen atoms with H 2 , and it dominates at A < 6000 A for cool 
white dwarfs, i.e. when H2 forms. Even though the Ly a opacity mainly supresses the flux 
in the ultraviolet, it causes a significant redistribution of flux toward longer wavelengths. 

We perform three separate sets of fits. The first set of fits involves pure hydrogen models 
with T eff and log g as free parameters. The second set of fits uses mixed H/He models and 
assumes that the helium-to-hydrogen ratio is a free parameter as well. Both sets of fits use 
models with T e g > 3500 K and log g = 7 — 9. F inally, the last set of fits were performed with 
blackbody SEDs. IKowalski fc Saumonl (120061 ) demonstrated that pure helium atmosphere 
white dwarfs have SEDs similar to blackbodies, therefore if any of our targets have pure 
helium atmospheres, they should be best explained by the third set of fits. 



4. Results 

Figure 1 displays the optical and infrared SEDs (error bars) for 44 cool white dwarfs with 
T eff < 6000 K. Three of these stars (plotted at the end of the figure) are from the SDSS, and 
their photometry is in the SDSS ugriz and 2MASS JHK S system. The remaining 41 stars 
have BVRIJHK photometry available. The expected fluxes from synthetic photometry 
of pure hydrogen and mixed H/He white dwarf model atmospheres are shown as open and 
filled circles, respectively. Since the pure hydrogen or mixed H/He models reproduce the 
SEDs better than the blackbody SEDs, the blackbody fits are not shown in any of the panels 
except for WD 2251—070. This star is a DZ white dwarf. Excluding the B-band flux which 
is affected by the presence of metals in the photosphere, WD 2251—070 is the only object 
for which a blackbody SED is a better fit to the observed SED than pure hydrogen or mixed 
H/He models. 
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The differences between the pure hydrogen and mixed H/He models are not statistically 
significant for the majority of our targets. For example, the predicted SEDs for the pure H 
and mixed H/He solutions for WD 1656—062 are almost identical. One cannot differentiate 
between these two models based on the SED alone. WD 1656—062, like 19 other WDs 



in our sample, is a DA white dwarf. iBergeron et al.l (120011 ) found that the pure hydrogen 
models are in good agreement with Ha spectroscopy of most stars in their (and our) sample. 
Therefore, we adopt the pure hydrogen atmosphere solutions for all DAs. 

A 5% uncertainty in the K— band flux corresponds to a helium detection limit of He/H = 
0.04, 0.36, and 2.14 for 4000, 5000, and 6000 K white dwarfs, respectively. Out of the 23 DC 
white dwarfs in our sample, 12 are best explained with mixed atmosphere models that have 
He/H ratios above this detection limit. Hence, we classify these white dwarfs as mixed H/He 
atmosphere objects. A few examples of the superior fits by mixed H/He fits include WD 
0029-032, WD 0747+073B, WD 1313-198, WD 1444-174, and SDSS J0753+4230. The 
differences between the best-fit pure hydrogen and mixed H/He models are not significant 
for the remaining 11 DCs, and we assign pure hydrogen composition to them. The physical 
parameters of the best fit atmosphere solutions as well as the spectral types of our targets 
are presented in Table 2. 

Only one star in Ta ble 2 is assigned a pur e heliu m composition. This result was expected 
considering the work of iKowalski fe Saumonl (120061 ). which has satisfactorily demonstrated 
that the optical and near-infrared colors of cool white dwarfs are explained fairly well with 
the new p ure hydrogen models. Yet th is result is surprising from an evolutionary point of 



view, since 



Tremblay fc Bergeron! (120081 ) demonstrated that 45% of white dwarfs in the range 



7000 K > T c ff > 5000 K have helium- rich atmospheres (with possible traces of hydrogen). 
The mixed H/He composition stars make up 27% of our sample. However, the helium to 
hydrogen ratio for these stars is on the order of unity. Our result presents an important 
challenge to understanding the spectral evolution of white dwarfs. 

In any case, the choice of pure hydrogen or mixed H/He models do not have an important 
effect on understanding the mid-infrared flux distributions of cool white dwarfs. Figure 1 
demonstrates that the main difference between the best-fit pure hydrogen and mixed H/He 
models is in the optical and near-infrared, and both sets of models agree in the mid-infrared. 
Moreover, it shows that our grid of models reproduce the observed SEDs well in the entire 
range from 0.4/xm to 8/zm. We do not find any significant mid-infrared flux deficits for these 
44 cool white dwarfs. We now check to se e how well the models are able to reproduce the 
observations. iTremblay fc Bergeron! (120071 ) found that the 12 cool white dwarfs that they 
analyzed were on average fainter by 2% in th e 4.5 Aim band, wh ereas they could fit the 8 
/im photometry fairly well with their models. iDebes et al.l (120071 ) and iFarihi et al.l (j2008al ) 
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found that the models reproduce the observed mid-infrare d fluxes within 5%, slightly larger 
than the IRAC calibration uncertainty of 3% measured by lReach et al.l (120051 ). 



Figure 2 presents the ratio of observed to predicted fluxes in all 11 filters from our 
analysis using the pure hydrogen atmosphere fits. The IRAC 3.6, 4.5, 5.8, and 8 /im bands 
are labeled as II, 12, 13, and 14, respectively. Figure 3 shows the same ratios for the best-fit 
models to the SEDs (31 pure hydrogen and 12 mixed H/He atmosphere solutions given in 
Table 2). A comparison of Figure 2 with Figure 3 shows that the use of mixed H/He models 
for the 12 DCs improves the agreement between the observations and models. Instead of 
looking at how well the models can reproduce the observations just in the mid-infrared, 
testing the models in all filters gives us a better idea of how well the models work in genera l. 
The photometric uncertainties in VRI are 3% and in BJHK are 5% (IBergeron et al.ll200ll ). 
Our photometric m ethods are consistent with those used for the IRAC absolute calibration 
(IReach et al.l 120051 ) . and therefore the mid-infrared fluxes should be accurate to 3%. 



Table 3 presents the mean differences between the observed photometry and models 
for all filters. The pure hydrogen models are consistent with the observations within 6%. 
Excluding the IRAC 8/i band which has the lowest signal-to- noise ratio photometry, the 
best fit models are consistent with observations within 5%. Given the 3-5% uncertainty in 
the observations, we conclude that the optical, near-infrared, and IRAC photometry of cool 
white dwarfs are consistent with white dwarf model atmospheres. 

Subtle systematic trends as a function of effective temperature are visible in Figure 2 
and 3. Even though individual points in these trends are not statistically off the one-to-one 
line, having a relatively large sample of stars enables us to see these subtle effects. The B— 
band and IRAC2 differences are perhaps the most convincing of these trends. Understanding 
these trends is useful for improving the current model atmosphere calculations. The observed 
trend in the B— band is especially important as the B— band flux is mostly suppressed by 
the Ly a red wing opacity in cool white dwarf atmospheres. Due to the uncertainty in the 
H2— H potential ene r gy su rface and transition dipole moment s, the Ly a profiles used in 
Kowalski fc Saumonl ( 20061 ) models are reliable up to 6000 A. iKowalski fc Saumonl ( 20061 ) 
extrapolated the Ly a profiles for longer wavelengths. The Ly a absorption dominates up to 
4500 A for T c g = 5000 K white dwarfs. However, it is important up to 6000 A and 7000 A 
for T e fj = 4500 K and 4000 K, respectively. Therefore, assuming that the other absorption 
mechanisms are correctly described in the models, the observed drift with T e g may indicate 
the need for extension of calculated Ly a profiles beyond A > 6000 A. However, we do not 
exclude other mechanisms that may cause the observed systematic trends. 



Figure 4 shows four white dwarfs with Spitzer photometry that are perhaps affected by 
nearby bright stars. Similar to the other 44 stars presented in previous figures, the SEDs 
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for these stars are also consistent with model predictions, and we do not find significant 
mid-infrared flux deficits. 

The optical and infrared photometry and our IRS spectrum of the brightest target in our 
sample, the DA white dwarf Wolf 489, are presented in Figure 5. The best fit model spectra 
with pure H and mixed H/He compositions are shown as solid and dashed lines, respectively. 
The main difference between these two models are in the optical and they predict the same 
flux distribution in the mid-infrared. The IRS spectrum is s hown as a jagged line with error 



bars. Since Wolf 489 has a proper motion of ~ 3.9" yr -1 ( iLepine &: Shard 120051 ) . the IRS 
centering was not perfect, and the measured fluxes from the spectrum were a factor of six 
lower than the IRAC photometry. The centering problem resulted in lower signal-to-noise 
ratio data and low level fluctuations are visible in the spectrum. Since the IRS spectrum 
covers the IRAC 8/xm band, we normalized it to match the 8/im photometry. A comparion 
of the best fit models and the observations shows that the observed IRS spectrum is entirely 
consistent with T e g ~ 5000 K white dwarf models over the 5-15 /zm range, further confirming 
our conclusion that the SEDs of cool white dwarfs are consistent with model predictions. 



4.1. Ultracool White Dwarfs 



White dwarfs with T e g thought to be below 4000 K are classifi ed as ultracool. Te n years 
after the discovery of the first ultraco ol white dwarf LHS 3250 ([Harris et al.lll999l ). there 



are now more than a dozen known (see Gates et al. 2004: Harris et al 



2008 



and references 

therein). However, only a few of them have been studied in de tail using optical and near - 
infrared photometry and trigonometric parallax measurements (IBergeron fc Leggettl 120021 ) . 
The main difference between cool and ultracool white dwarfs is that several of these ultracool 
white dwarfs show significant flux deficits in the near-infrared and even in the optical, which is 
most probably caused by the strong CIA absorption from H2 in a dense, helium dominated 
atmosphere. The SEDs of these white dwarfs cannot be fit with any model atmospheres 
currently available. 

LHS 1126 is not classi fied as an ultracool w hite dwarf since its optical colors are con- 
sistent with T e fr = 5400 K (IBergeron et al.lll994j ). However, it also displays significant flux 
deficits in the near- and mid-infrared. A remaining question is whether all white dwarfs that 
display significant flux deficits in the near-infrared also display deficits in the mid-infrared 
or not. The importance of this question is that current CIA opacity calculations do not pre- 
dict absorption in the mid-infrared. Therefore, finding several ultracool white dwarfs with 
mid-infrared flux deficits would mean that there is a missing opacity source in current white 
dwarf atmosphere models. 
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Three ultracool white dwarfs were included in our sample, LHS 3250, CE 51, and WD 
0346+246. Unfortunately, WD 0346+246 was blended with a brighter source in the IRAC 
images and we could not reco ver its photometry. Spitz er photometry of CE 51 is also affected 



by a very bright nearby star. iRuiz &: Bergeron! (120011 ) assigned a temperature of 2730 K to 



CE 51 based on its optical photometry. Extending our models down to 2500 K, and including 
the IRAC photometry in our fits does not change the best-fit temperature much. We find 
a best-fit solution of 2660 K for a pure hydrogen composition. However, the best-fit model 
is not a good representation of the optical and mid-infrared photometry. Since our IRAC 
photometry is questionable for this star, we cannot draw any other conclusions at this time. 

LHS 3250 was close to an extended source in our IRAC images, however we are able 
to recover the photometry using a 2 pixel aperture. Figure 6 shows LHS 3250 images from 
1988 to 2007. Our Spitzer IRAC images were taken in August 2006, 9 months before the 
/—band image from the MDM 2.4m telescope was taken. LHS 3250 has a proper motion 
of cosS,fis) = (—0.541, +0.148 mas yr -1 ; Munn et al. 2004) , therefore the position 
measured from the J— band images should be essentially the same as the position in our 
IRAC images. This position is marked with open circles in the MDM and IRAC images; 
LHS 3250 is clearly detected. 

LHS 3250 is quite faint in the mid-infrared. From the BCD frames, we measure average 
fluxes of 7.0 ± 1.7, 5.8 ± 2.3, 8.1 ± 4.9, and 13 ± 6.4 /iJy in the 3.6, 4.5, 5.8, and 8 fim 
bands, respectively. Using the combined mosaic images, we measure 7.0, 6.0, 6.1, and 10.7 
//Jy in the same bands, respectively. Photometry from the BCD frames and the mosaic 
frames are consistent within the errors. The LHS 3250 SED compared to the best fit mixed 
H/He atmosphere model is shown in the top panel in Figure 7. The effect of CIA opacity 
is clearly seen in the models; there is a big dip in the predicted fluxes up to about 3 //m 
due to H 2 1-0 dipole absorption centered at 2.4 /im. The models obviousl y fail to reproduce 



the observations in all ban ds. As in previous attempts by other authors ([Harris et al.lll999 



Bergeron &: Leggettl 120021 ). the SED cannot be explained with current white dwarf model 



atmospheres. 

The bottom panel in Figure 7 compares the LHS 3250 SED with LHS 1126. We normal- 
ized the LHS 1126 SED to match LHS 3250 in the B— band. Since LHS 1126 is more than 
an order of magnitude brighter, its infrared fluxes are more accurate. The infrared portion 
of the LHS 1126 SED is best fit with a power law index of —1.99. On the other hand, the 
infrared portion of the LHS 3250 SED, excluding the 5.8 and 8 /im photometry due to large 
error bars, fits a power law index of —2.07. However, the power law fit is not as good a fit 
as for LHS 1126. The fluxes measured between 3.6 and 8 /im are consistent with each other 
within the errors. 
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5. Discussion 
5.1. Cool White Dwarf Atmospheres in the Infrared 

We now have extend ed the earlier work on Spitzer observations of cool white dwarfs 



( IKilic et al.l l2006al . l2008al ) down to 4000 K using a well studied sample of nearby white 
dwarfs. Our fits to the BVRIJHK and 3.6 — 8/mi photometry of ~50 cool white dwarfs 
show that current state of the art white dwarf model atmospheres are able to reproduce the 
observations within 5%. The only exception to this is LHS 1126, which shows significant flux 
deficits in the mid-infrared compared to the models. In addition, our models are not able 
to explain the observed optical and infrared SED of the ultracool white dwarf LHS 3250. 
More theoretical work on understanding the properties of dense helium- rich, mixed H/He 
white dwarf atmospheres is probably needed to understand the spectra of stars like this. A 
comparison of the SEDs for LHS 1126 and LHS 3250 shows that a power law may explain 
the infrared flux distributions of both stars, however, the LHS 3250 SED shows evidence of 
a flat or increasing mid-infrared flux distribution, and therefore its SED remains puzzling. 
Obtaining higher quality mid-infrared data for LHS 3250 and extending the wavelength 
coverage further to the red will be important in understanding LHS 3250 and other ultracool 
white dwarf atmospheres. This will be possible with the Mid-Infrared Instrument (MIRI) 
on the James Webb Space Telescope. 



5.2. Infrared Excess from Dust Disks 



Previous ground- and space-based studies of metal-rich white dwarfs showed that ~20% 
of hydrog en atmosphere whi te dwarfs with T e g > 7400 K display infrared excesses from dust 



disks (see IKilic et al.l l2008bl . and references therein). However, the majority of these disks 



are found around stars hotter than 10000 K and younger than 1 Gyr. The fraction of dusty 
white dwarfs goes down to 1-3% for previous Spitzer o bservations of ^ 200 white dwarfs with 



somewhat less metallicity bias (IFarihi et al.l l2008bl ) . iMullally et al.l (120071 ) found 2 dusty 



stars out of 124 white dwarfs, corresponding to a fraction of l^]^ A ssuming optically 



thick and flat disk mo dels, these disks are within IRq of white dwarfs (jJura et al.l 120071 ; 
von Hippel et al.l 120071 ) . An important step in understanding the formation and evolution 
of these disks is the disk lifetimes and whether the dust can survive for the lifetime of the 
white dwarf or not. 

Previous Spitzer observations of white dwarfs targeted mostly stars with T e g > 6000 K 
and ages younger than 2 Gyr. Extending the temperature limit down to 4000 K correspond s 
to extending the ages up to 9 Gyr for typical 0.6 M & white dwarfs ( iBergeron et al.lll995l ). 
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Excluding the ultracool white dwarfs from our sample, none of the 48 cool white dwarfs 
(Figure 1 and 4) in our sample shows excess infrared radiation from disks. However, cooler 
stellar effective temperatures would mean that if disks are present, they will be substantially 
cooler than the disks around warmer white dwarfs. For example, if we take the disk around 
G29-38 as a prototype and put it around a 5000 K white dwarf, the disk would have an inner 
temperature of ~ 600 K and an outer temperature of 350 K. 



Figure 8 shows several flat disk models (IJural 120031 ) for a 5000 K white dwarf at 20 
pc. The solid line shows the total flux expected from the photosphere of a 5000 K white 
dwarf and a disk with an inclination angle 45° and the same inner-outer radii as G29-38. 
The short-dashed line shows the flux expected from a disk further away from the star, with 
parameters similar to G 166— 58 (IFarihi et al.ll2008al ). The other dashed and dashed-dotted 
lines show the differences in disk emission depending on the inner radii of the disks. Here 
we only explo re the disks out to 1.2 as it correspon ds to the Roche limit for typical 
white dwarfs (jvon Hippel et al.l 120071 ; IFarihi et al.l l2008al ) . Adopting an accuracy of 5% in 
the Spitzer IRAC bands, the solid horizontal line shows 3a excess above the photospheric 
emission at 8//m. It is evident from this figure that, if disks survive for several billion years 
or if they are replenished, they should be detected in our IRAC observations if the inner 
radii of the disks are smaller than 0.5 R@. Colder disks with larger inner radii would be 
visible at longer wavelengths. 

Non-detection of disks around our 48 targets means that the frequency of warm disks 
around cool white dwarfs is 0.8]t 'g %, still consistent with the frequency of disks for warmer 
white dwarfs. 



MK thanks D. Winget and S. Kenyon for useful discussions. Support for this work was 
provided by NASA through the Spitzer Space Telescope Fellowship Program, under an award 
from Caltech, and also award projects NBR 1359963 and 1288635 issued by JPL/Caltech to 
Ohio State University and the University of Texas. PMK thanks D. Marx and the theoretical 
chemistry group at Ruhr-Universitat Bochum for their hospitality and support. This work 
is based in part on observations made with the Spitzer Space Telescope, which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology under NASA contract 
1407. 



REFERENCES 

Bergeron, P., Ruiz, M.-T., Leggett, S. K., Saumon, D., & Wesemael, F. 1994, ApJ, 423, 456 



Bergeron, P., Wesemael, F., & Beauchamp, A. 1995, PASP, 107, 1047 
Bergeron, P., Ruiz, M. T., & Leggett, S. K. 1997, ApJS, 108, 339 
Bergeron, P., Leggett, S. K. , & Ruiz, M. T. 2001, ApJS, 133, 413 
Bergeron, P., & Leggett, S. K. 2002, ApJ, 580, 1070 

Bergeron, P., Ruiz, M. T., Hamuy, M., Leggett, S. K., Currie, M. J., Lajoie, C.-P., & Dufour, 
P. 2005, ApJ, 625, 838 

Cutri, R. M et al. 2003, VizieR Online Data Catalog, 2246 

Debes, J. H., Sigurdsson, S., & Hansen, B. 2007, AJ, 134, 1662 

Farihi, J., Zuckerman, B., & Becklin, E. E. 2008a, ApJ, 674, 431 

Farihi, J., Jura, M., & Zuckerman, B. 2008b, ApJ, submitted 

Fazio, G. G. et al. 2004, ApJS, 154, 10 

Gates, E. 2004, ApJ, 612, 129 

Hall, P. B., & Maxwell, A. J. 2008, ApJ, 678, 1292 

Hambly, N. C, Smartt, S. J., & Hodgkin, S. T. 1997, ApJ, 489, L157 

Hansen, B. M. S. 1998, Nature, 394, 860 

Hansen, B. M. S., Brewer, J., Fahlman, G. G., Gibson, B. K., Ibata, R., Limongi, M., Rich, 
R. M., Richer, H. B., Shara, M. M., & Stetson, P. B. 2002, ApJ, 574, L155 

Hansen, B. M. S., et al. 2007, ApJ, 671, 380 

Harris, H. C, Dahn, C. C, Vrba, F. J., Henden, A. A., Liebert, J., Schmidt, G. D., & Reid, 
I. N. 1999, ApJ, 524, 1000 

Harris, H. C, et al. 2006, AJ, 131, 571 

Harris, H. C, et al. 2008, ApJ, 679, 697 

Higdon et al, PASP 116, 975-984, 2004 

Houck, J. R., et al. 2004, ApJS, 154, 18 

Jura, M. 2003, ApJ, 584, L91 



Jura, M., Farihi, J., Zuckerman, B., k Becklin, E. E. 2007, A J, 133, 1927 

Kilic, M., von Hippel, T., Mullally, F., Reach, W. T., Kuchner, M. J., Winget, D. E., k 
Burrows, A. 2006, ApJ, 642, 1051 

Kilic, M., et al. 2006, AJ, 131, 582 

Kilic, M., Kowalski, P. M., Mullally, F., Reach, W. T., k von Hippel, T. 2008, ApJ, 678, 
1298 

Kilic, M., Farihi, J., Nitta, A., k Leggett, S. K. 2008, AJ, 136, 111 

Koester, D. k Wolff, b. 2000, A&A, 357, 587 

Kowalski, P. M. k Saumon, D. 2004, ApJ, 607, 970 

Kowalski, P. M. 2006, ApJ, 641,488 

Kowalski, P. M., & Saumon, D. 2006, ApJ, 651, L137 

Kowalski, P. et al 2007, Phys. Rev. B. 2007, 76 

Leggett, S. K., Ruiz, M. T., k Bergeron, P. 1998, ApJ, 497, 294 

Lepine, S., k Shara, M. M. 2005, AJ, 129, 1483 

Mestel, L. 1952, MNRAS, 112, 583 

Mullally, F., Kilic, M., Reach, W. T., Kuchner, M. J., von Hippel, T., Burrows, A., k 
Winget, D. E. 2007, ApJS, 171, 206 

Munn, J. A., et al. 2004, A J, 127, 3034 

Reach, W. T., et al. 2005, PASP, 117, 978 

Ruiz, M. T., k Bergeron, P. 2001, ApJ, 558, 761 

Saumon, D. k Jacobson, S. B. 1999, ApJ, 511, 107 

Tremblay, P.-E., k Bergeron, P. 2007, ApJ, 657, 1013 

Tremblay, P.-E., k Bergeron, P. 2008, ApJ, 672, 1144 

von Hippel, T., Kuchner, M. J., Kilic, M., Mullally, F., k Reach, W. T. 2007, ApJ, 662, 544 



-14- 



Winget, D. E., Hansen, C. J., Liebert, J., van Horn, H. M., Fontaine, G., Nather, R. E., 
Kepler, S. O., & Lamb, D. Q. 1987, ApJ, 315, L77 



This preprint was prepared with the AAS iATgX macros v5.2. 



- 15 - 



Table 1. IRAC Photometry of Cool White Dwarfs 



Object 


Name 




3.6/im 
(mJy) 






4.5^m 
(mJy) 






5.8^m 
(mJy) 






8%m 
(mJy) 




WU UU11 — 399 


TAA1 A OAOT 

J0014— 3937 


0. 


1186 


± 


0. 


AAO A 

.0039 





.0798 


± 





AAOA 

.0029 


A 




.0597 


± 


0. 


A ATO 

.0073 


0. 


.0330 ± 





AI O O 

.0132 


WU 0029— 032 


t tto 1 nno 

LHS1093 


A 

0. 


.2591 


± 


0, 


AAO £ 

0085 


A 




.1824 


± 





A A/ 1 T 

.0067 


A 

0. 


.1274 


± 


0. 


AOA A 

.(Jzd4 


A 

0. 


.0502 ± 





AI AO 

.0193 


~\htt\ nn^o oot 

WU 0042— 337 


TT'O E 1 E A 

b 351—50 


A 

0. 


0613 


± 





A AOO 


A 




,0426 


± 





A A 1 O 

.0018 


A 

0. 


.0283 


± 


0. 


A 1 OA 

0120 


A 

0. 


.0216 ± 





AAOA 

,0039 


~\HTT\ AA /i O OAT 

WD 0048—207 


T TTC1 1 C O 

LHS1158 


A 

0. 


.2308 


± 


0. 


AAO 1 

0081 


A 




.1550 


± 





A A C T 

.0057 


A 

0. 


,0965 


± 


0. 


AAO O 

.0088 


A 

0. 


.0453 ± 





A A 1 O 

.0018 


WD 0102+210.2 


T TTOrAO/1 

LHS5024 


0. 


.1129 


± 


0. 


.0036 





.0779 


± 





.0028 


0. 


.0505 


± 


0. 


.0103 


0. 


.0249 ± 





,0241 


WD 0102+210.1 


T TTCTAAO 

LHS5023 


0. 


.1018 


± 


0, 


0033 





.0715 


± 





.0026 


0. 


,0536 


± 


0. 


.0104 


0. 


.0224 ± 





.0161 


~\T7"T~* AI 1 n A 1 O 

WD 0112—018 


T TT0 1 Ol A 

LHS1219 


0. 


.1383 


± 





.0044 





.0908 


± 





.0032 


0. 


.0598 


± 


0. 


.0157 





.0324 ± 





,0144 


WD 1)117—145 


T TTC1 OO O 

LHS1233 


A 

0. 


.3033 


± 





AAAO 

.0098 


A 




.2032 


± 





A A T 1 

.0071 


A 




.1239 


± 


0. 


.0223 


A 




.1083 ± 





AOA A 

,0266 


"\ ~x 7~T\ a 1 1 i /i a 1 

WD 0121+401 


G133— 8 


0. 


.2216 


± 


0, 


.0074 





.1465 


± 





.0055 


0. 


.0831 


± 





.0098 


0. 


.0567 ± 





.0102 


~\ I 7"T~* A1 OO f A O 

WD 0133—548 


"ITS 7TO A 1 r r - /I A 

WD0135— 546 


0. 


.1044 


± 


0, 


.0034 





.0686 


± 





.0024 


0. 


.0457 


± 


0. 


.0050 


0. 


.0263 ± 





.0049 


~\TS7T~\ AO AO Arr 

WD 0202—055 


7TO AO A E A E O 


A 

0. 


.1057 


± 


0, 


AAO tZ 

.0035 


A 




.0732 


± 





AAO O 

.0033 





.0399 


± 


0. 


.0078 





.0277 ± 





,0066 


~\*7T~Y AOOO 1 A/IOO. 

WD 0222-\-bA& 


T TTO 1 /I AC 

LHbl405 


A 

0. 


.1358 


± 


0, 


ao c: 

0352 


A 




.0951 


± 





AO A A 

.(Jz44 


A 

0. 


.0653 


± 


0. 


A 1 tr A 

.0159 


A 

0. 


.0401 ± 





A 1 A A 

.0100 


~\XTT\ A O O A 1/1/1 

WD Uz3U— 144 


T TTO 1 A 1 E 

LHbl415 


A 

0. 


.7045 


± 





AOO A 

.0220 


A 




.4654 


± 





A 1 C O 

.0152 


A 




.3185 


± 


0. 


A 1 A A 

Ulo4 


A 




.1754 ± 





A 1 OA 

,0180 


~\*7T~\ AOtTT 1 A0 1 

WD 0357+081 


T TTO 1 A 1 T 

LHS1617 


A 

0. 


.6134 


± 


0. 


A 1 AO 

.0193 


A 




.4106 


± 





A 1 Of 

.0136 


A 




.2526 


± 





A 1 tr A 

.0159 


A 

0. 


.1416 ± 





A 1 A A 

.0140 


T T * 1 \ A ,1 AT 1 1 A T 

WD 0407+197 


T TTO 1 AO/ 1 

LHS1636 


0. 


,1354 


± 


0, 


.0050 





.0815 


± 





.0040 





.0504 


± 


0. 


.0244 





.0376 ± 





.0289 


T ~X TT~\ A /I O O 1 A A A 

WD 0423+044 


T TTO 1 O T A 

LHS1670 


0. 


.2655 


± 


0. 


.0088 





.1789 


± 





.0069 


0. 


.1185 


± 





.0123 


0. 


.0785 ± 





.0189 


~\T7"T~* AT AO 1 T /I 

WD 0503—174 


T TTO 1 TO /I 

LHS1734 


0. 


,5840 


± 


0. 


.0187 





.3797 


± 





.0125 


0. 


.2536 


± 


0. 


.0164 





.1455 ± 





.0142 


~\*7T~Y ACTCT1 i AfOCL 

WD 1)551+408 


T TTO 1 OA 1 

LHS1801 


A 

0. 


.2527 


± 


0. 


AAO A 

.0089 


A 




.1608 


± 





.0068 


A 




.0979 


± 


0. 


A 1 O A 

.0134 


A 




.0510 ± 





AO A C 

,0205 


~\ T ' I \ A/" 1 C T 1 O OA 

WD 0657+320 


T TTO 1 O OA 

LHS1889 


0. 


,4238 


± 


0. 


.0135 





.2861 


± 





.0097 





.1969 


± 


0. 


.0140 


0. 


.1279 ± 





.0169 


WD 0727+482AB 


/~1 1 AT TA 

G107— 70 


2. 


,5598 


± 





.0782 


1 


.7041 


± 





.0525 


1 


.1482 


± 


0. 


.0395 





.6676 ± 





,0256 


~\ ~X TT~\ AT 1 O O /I A 

WD 0743—340 


vB3 






± 









.3585 


± 





.0125 






± 






0. 


.1044 ± 





.0080 


Tim A7/I7 1 A70 A 

WD 0747+073A 


T TTO OO A 

LHS239 


A 

0. 


.3987 


± 


0. 


A 1 OT 

0127 


A 




.2739 


± 





AAA A 

.0094 


A 

0. 


.1668 


± 


0. 


A 1 A 1 

.1)141 


A 

0. 


.1017 ± 





AO AO 

.0202 


TX7T\ AT /IT 1 A TO TO 

WD 0747+073B 


T TTOO A (~\ 

LHS240 


A 

0. 


.4552 


± 


0. 


.0144 


A 




.3085 


± 





A 1 AO 

.0103 


A 




.1990 


± 


0. 


A1 AO 

.0168 


A 




.1048 ± 





Ai rr 

.0155 


SDSS J0753+4230 


J0753+4230 


0. 


,2028 


± 


0. 


.0068 





.1427 


± 





.0056 


0. 


0963 


± 


0. 


.0093 





.0565 ± 





.0189 


WD 0851-246" 


CE51 


0. 


.1010 


± 


0. 


.0066 





.0790 


± 





.0064 





.0677 


± 


0. 


.0164 




... ± 






WD 1022+009 


LHS282 


0. 


,0875 


± 


0. 


0029 





.0557 


± 





.0023 


0, 


,0274 


± 





.0093 


0. 


.0183 ± 





.0126 


WD 1108+207 


LHS2364 


0. 


.1892 


± 





.0067 





.1197 


± 





.0055 





,0955 


± 


0. 


0133 





.0310 ± 





.0346 


WD 1136-286 


ES0439-26 


0. 


.0144 


± 


0. 


.0010 





.0077 


± 





.0014 





,0186 


± 


0. 


.0110 





.0031 ± 





.0003 


WD 1153+135 


LHS2478 


0. 


.1631 


± 


0. 


0056 





.1050 


± 





.0046 


0, 


,0800 


± 





.0200 


0. 


.0559 ± 





.0188 


WD 1247+550° 


LP131-66 


0. 


2353 


± 





.0082 





.1569 


± 





.0060 





.1025 


± 


0. 


.0111 





.0648 ± 





.0106 


WD 1257+037 


LHS2661 


0. 


.5908 


± 


0. 


.0187 





.3859 


± 





.0126 





.2668 


± 


0. 


.0163 


0. 


.1340 ± 





.0161 


WD 1300+263 


LHS2673 


0. 


.0749 


± 


0. 


,0026 





.0517 


± 





.0023 


0, 


.0399 


± 


0. 


.0061 


0. 


.0121 ± 





.0103 


WD 1313-198 


LHS2710 


0. 


2037 


± 


0. 


,0069 





.1336 


± 





.0052 





.1003 


± 


0. 


.0147 





.0722 ± 





.0240 


SDSS J1313+0226 


J1313+0226 


0, 


,1327 


± 


0. 


,0044 





.0901 


± 





.0031 


0. 


.0615 


± 





.0069 


0. 


.0379 ± 





.0128 


WD 1334+039 


Wolf489 


2. 


,6091 


± 


0. 


,0805 


1 


.7405 


± 





.0548 


1 


.1851 


± 


0. 


.0477 





.6677 ± 





.0329 


WD 1345+238 


LP380-5 


1. 


,1767 


± 


0. 


0363 





.7908 


± 





.0248 


0. 


.5380 


± 





.0226 


0. 


.3165 ± 





.0185 



Table 1 — Continued 



Object 


Name 


3.6fim 




5.8/im 


8.0^m 






(mJy) 


(mJy) 


(mJy) 


(mJy) 



WD 1346+121 


LHS2808 





.1001 


± 





.0035 


0. 


,0654 


± 


0, 


.0043 


0. 


,0411 


± 





.0125 





.0330 


± 


0. 


,0036 


WD 1444-174 


LHS378 





.4579 


± 





.0146 


0. 


,3142 


± 


0, 


.0106 


0, 


,2212 


± 





.0174 





.1589 


± 


0. 


,0162 


WD 1602+010 


LHS3151 





.1701 


± 





.0058 


0. 


,1151 


± 


0, 


.0047 


0, 


,0671 


± 





.0114 





.0498 


± 


0. 


.0164 


WD 1653+630 


LHS3250 





.0070 


± 





.0017 


0. 


,0058 


± 


0. 


,0023 


0. 


.0081 


± 





.0049 





.0130 


± 


0. 


.0064 


WD 1656-062 


LP686-32 





.1826 


± 





.0069 


0. 


.1203 


± 


0, 


.0055 


0, 


.0746 


± 





.0186 





.0203 


± 


0. 


.0156 


WD 1820+609 


G227-28 


1 


.0795 


± 





.0333 


0. 


.7283 


± 





,0229 


0, 


,4773 


± 





.0198 





.2700 


± 


0. 


.0201 


WD 2002-110 


LHS483 





.3407 


± 





.0115 


0. 


.2390 


± 


0. 


.0090 


0, 


.1359 


± 





.0157 





.0956 


± 


0. 


.0205 


WD 2048+263 


G187-8 





.9917 


± 





.0307 


0. 


.6574 


± 


0. 


.0208 


0, 


.4573 


± 





.0206 





.2591 


± 


0. 


.0146 


WD 2054-050" 


vBll 





.5335 


± 





.0235 


0. 


.3487 


± 


0. 


.0270 


0. 


.1060 


± 





.0820 





.0778 


± 


0. 


0266 


SDSS J2116-0724 


J2116-0724 





.1442 


± 





.0051 


0. 


.0991 


± 





.0042 


0. 


.0735 


± 





.0327 





.0532 


± 


0. 


.0214 


WD 2248+293 


G128-7 





.7773 


± 





.0242 


0, 


.5090 


± 


0, 


,0165 


0. 


.3545 


± 





.0184 





.1781 


± 


0. 


.0142 


WD 2251-070 


LP701-29 


1 


.3402 


± 





.0412 


0. 


.9103 


± 


0, 


,0287 


0. 


.6108 


± 





.0273 





.3791 


± 


0. 


.0252 


WD 2316-064 


LHS542 





.1237 


± 





.0040 


0. 


.0799 


± 


0. 


,0030 


0. 


,0596 


± 





.0085 





.0283 


± 


0. 


.0165 


WD 2343-481 


WD2346-478 





.1225 


± 





.0039 


0. 


.0752 


± 


0, 


.0028 


0. 


.0452 


± 





.0052 





.0323 


± 


0. 


.0116 


WD 2345-447 


ES0292-43 





.0948 


± 





.0031 


0. 


.0627 


± 


0, 


,0023 


0. 


,0431 


± 





.0067 





.0319 


± 


0. 


.0121 



Note. — (a)- The photometry is likely to be affected by nearby bright star, (b)- WD 0727+482A and WD 
0727+482B are unresolved in IRAC images. The photometry given here is the combined flux from both stars. 
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Table 2. Atmospheric Parameters of Cool White Dwarfs 



Object 


Spectral 
Type 


Comp 


Tcff 
(K) 


log g 


Hc/H 


WD 0011-399 


DC 


H 


4700 


8 




WD 0029-032 


DC 


Hc/H 


4590 


7.96 


0.3 


WD 0042-337 


DC 


H 


4040 


8 




WD 0048-207 


DA 


H 


5190 


8 




WD 0102+210.2 


DC 


H 


4910 


8 




WD 0102+210.1 


DA 


H 


5320 


8 




WD 0112-018 


DA 


H 


5410 


8 




WD 0117-145 


DA 


H 


5070 


7.75 




WD 0121+401 


DA 


H 


5330 


7.94 




WD 0133-548 


DC 


Hc/H 


4660 


8 


0.2 


WD 0202-055 


DC 


H 


4150 


8 




WD 0230-144 


DA 


H 


5470 


8.15 




WD 0357+081 


DA 


H 


5490 


8.07 




WD 0407+197 


DC 


H 


5250 


8 




WD 0423+044 


DA 


H 


4950 


8 




WD 0503-174 


DA 


H 


5410 


7.72 




WD 0657+320 


DA 


H 


4930 


8.09 




WD 0743-340 


DC 


H 


4630 


8.15 




WD 0747+073A 


DA 


H 


4500 


7.97 




WD 0747+073B 


DC 


Hc/H 


4700 


7.97 


0.8 


SDSS J0753+4230 


DC 


He/H 


4540 


8 


0.7 


WD 1022+009 


DA 


H 


5400 


7.65 




WD 1108+207 


DC 


H 


4760 


8.18 




WD 1136-286 


DC 


H 


4630 


9.12 




WD 1153+135 


DC 


Hc/H 


4830 


8 


0.3 


WD 1257+037 


DA 


H 


5640 


8.22 




WD 1300+263 


DC 


H 


4360 


8.16 




WD 1313-198 


DC 


Hc/H 


5110 


8.24 


1.7 


SDSS J1313+0226 


DC 


H 


4230 


8 




WD 1334+039 


DA 


H 


5030 


8.01 




WD 1345+238 


DA 


H 


4690 


7.91 




WD 1346+121 


DC 


He/H 


4710 


8 


0.9 


WD 1444-174 


DC 


Hc/H 


4820 


8.31 


0.7 


WD 1602+010 


DC 


H 


4880 


8 




WD 1656-062 


DA 


H 


5550 


8.09 




WD 1820+609 


DA 


H 


4900 


7.97 
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Table 2 — Continued 



Object 


Spectral 


Comp 




log g 


Hc/H 




Type 




(K) 






WD 2002-110 


DC 


Hc/H 


4640 


8.23 


0.3 


WD 2048+263 


DA 


H 


5110 


7.65 




SDSS J2116-0724 


DC 


Hc/H 


4650 


8 


0.3 


WD 2248+293 


DA 


H 


5650 


7.65 




WD 2251-070 


DZ 


He 


4790 


Blackbody 




WD 2316-064 


DC 


H 


4700 


8.18 




WD 2343-481 


DA 


H 


5130 


8 




WD 2345-447 


DC 


Hc/H 


5330 


8.70 


2.0 
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Table 3. Mean % Differences Between the Observed Photometry and Models 





pure H 


pure H 


Best-fit 


Best-fit 


Filter 


Mean 


Wt-Mean 


Mean 


Wt-Mean 


B 


-5.7 ±3.7 


-6.0 ±0.7 


-4.9 ±3.7 


-5.2 + 0.7 


V 


+1.5 ±2.7 


+ 1.3 + 0.4 


+0.9 + 2.2 


+0.8 + 0.4 


R 


+3.8 ±3.2 


+3.6 + 0.5 


+2.8 + 2.8 


+2.7 + 0.5 


I 


+2.9 + 4.3 


+2.5 + 0.5 


+1.6 + 3.6 


+ 1.4 + 0.4 


J 


+6.0 + 4.4 


+5.7 + 0.8 


+5.0 + 5.3 


+4.5 + 0.8 


H 


-0.9 + 6.5 


-1.8 + 0.8 


+0.7 + 5.0 


+0.2 + 0.8 


K 


-0.9 + 6.6 


-1.9 + 0.8 


+1.8 + 5.1 


+ 1.4 + 0.8 


IRACf 


-2.0 + 4.2 


-2.5 + 0.5 


-1.1 + 3.6 


-1.5 + 0.5 


IRAC2 


-5.0 + 4.6 


-4.6 + 0.5 


-4.5 + 4.4 


-4.2 + 0.5 


IRAC3 


-O.f ±30.6 


-3.1 ± 1.3 


-0.1 + 30.6 


-3.0 ± 1.3 


IRAC4 


-5.7 ±21.1 


-11.4 + 1.6 


-6.3 + 20.3 


-11.6 ± 1.6 
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Fig. 1. — Spectral energy distributions of cool white dwarfs observed in our program. The 
observed fluxes are shown as error bars, whereas the expected flux distributions from pure 
hydrogen and mixed H/He atmosphere models are shown as open and filled circles, respec- 
tively. Instead of the mixed H/He atmosphere models, only the panel for WD 2251—070 
shows the predicted flux distribution for a blackbody SED as filled circles. 
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Fig. 1. — contd. 
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Fig. 2.- 



Ratio of the observed vs. model fluxes for the pure hydrogen models. 
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Fig. 3.- 



Ratio of the observed vs. model fluxes for the best-fit models. 
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Fig. 4. — SEDs of stars that have nearby bright stars in Spitzer images. 
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5 10 15 



Wavelength (fxm) 

Fig. 5. — The spectral energy distribution of Wolf 489 including the IRS spectrum from 5 
to 15 jum. The best-fit pure hydrogen and mixed H/He atmosphere models are shown as 
solid and dashed lines, respectively. The IRS spectrum is normalized to match the IRAC 
photometry in the 8 fim band. 
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Fig. 6.— LHS 3250 in June 1988 POSSII blue plates with North up and East to the left (top 
left panel), May 2007 MDM 2.4m telescope /—band image (top right panel), and August 
2006 Spitzer IRAC Channel 1 and 2 images (bottom left and right panels). 
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Fig. 7. — The LHS 3250 SED (error bars) compared to the best-fit mixed H/He model 
atmosphere (filled circles, top panel) and to LHS 1126 (open circles, bottom panel). The 
LHS 1126 SED is normalized to match LHS 3250 in the B— band. Power law fits to the 
infrared portion of the SEDs are shown as dashed lines. The power law fit to the LHS 3250 
SED does not include 5.8 and 8 jum photometry due to large error bars, and it has an index 
of -2.07. 
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Fig. 8. — The predicted flux distributions of a 5000 K white dwarf (dotted line) and hyp- 
thetical disks at different distances from the star. The inclination angle is fixed at 45°. The 
inner and outer radii of the disks are labeled in the figure. The thick solid line marks 3c 
detection limit at 8/im assuming a photometric accuracy of 5%. 



